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In Randall- Sundrum II (RS-II) braneworld model, it has been conjectured according to the 
AdS/CFT correspondence that brane- localized black hole (BH) larger than the bulk AdS curvature 
scale I cannot be static, and it is dual to a four dimensional BH emitting the Hawking radiation 
through some quantum fields. In this scenario, the number of the quantum field species is so large 
that this radiation changes the orbital evolution of a BH binary. We derived the correction to the 
gravitational waveform phase due to this effect and estimated the upper bounds on I by performing 
Fisher analyses. We found that DECIGO/BBO can put a stronger constraint than the current 
table-top result by detecting gravitational waves from small mass BH/BH and BH/neutron star 
(NS) binaries. Furthermore, DECIGO/BBO is expected to detect 10'' BH/NS binaries per year. 
Taking this advantage, we found that DECIGO/BBO can actually measure I down to ^ = 0.33/im 
for 5 year observation if we know that binaries are circular a priori. This is about 40 times smaller 
than the upper bound obtained from the table-top experiment. On the other hand, when we take 
eccentricities into binary parameters, the detection limit weakens to £ = 1.5/im due to strong de- 
generacies between I and eccentricities. We also derived the upper bound on I from the expected 
detection number of extreme mass ratio inspirals (EMRIs) with LISA and BH/NS binaries with 
DECIGO/BBO, extending the discussion made recently by Mc Williams [l|. We found that these 
less robust constraints are weaker than the ones from phase differences. 

PACS numbers: Valid PACS appear here 

I. INTRODUCTION 

Super-string theory suggests that our universe has more than four dimensions [H with extra dimensions being 
compactified in some way. One very famous and simple way is the Kaluza-Klcin compactification. The size of extra 
dimension t in this case is strongly bounded from particle physics experiments as £ < 10~^^cm. A new possibility 
opened after Arkani-Hamed et al. proposed a braneworld model (the ADD model) 0, Q (see also Ref. Q for the 
string realization of low scale gravity and braneworld models). They embedded a tension- less brane (on which we 
live) in a flat and compact bulk spacetime and assumed ordinary matters to be localized on it. In this case, gravitons 
are the only components that can propagate through the bulk. Since the law of gravity has been constrained only 
weakly by experiments, the size of extra dimensions can be relatively large in this model. Moreover, the ADD model 
can give an alternative way to explain the hierarchy problem between the Planck scale and the electroweak scale if 
the spacetime dimension is 6 and the size of compactified bulk is 1mm. 

A different type of braneworld models have been proposed by Randall and Sundrum H, 0] ■ In their first model 
(RS-I model) [6| , a positive tension brane and a negative tension brane (on which we live) give boundaries of a five 
dimensional bulk space with a negative cosmological constant. The unperturbed bulk is anti-de Sitter space whose 
metric is given by 

ds^ = e-^y'^rif.^dx^'dx'' + (1) 

where 77^^ represents the Minkowski metric with /i, v indices running from to 3, are the coordinates on the brane 
and y is the coordinate of the extra dimension, i is the AdS curvature length which characterizes the size of extra 
dimension. The location of each unperturbed brane is specified by a ?/ ^constant surface, and Z2 symmetry across 
the brane is assumed. This model can also solve the hierarchy problem by tuning the separation between two branes 
to be - 37£. 

In their second model (RS-II model) 0, it is assumed that we live on the positive tension brane. In this case, the 
position of the negative tension brane becomes almost irrelevant. Even if we send it to y = cx), the extra dimension 
is still effectively compactified thanks to the exponential warp factor e~'^yl^ in the metric. This model does not 
give any clue to the hierarchy problem, but it has a fascinating property: Although the model has a non-compact 
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extra dimension, four dimensional general relativity is approximately reproduced on the positive tension brane. The 
gravitational potential between two masses m and M with separation r on the brane becomes Q 



with G representing the effective gravitational constant. For the power-law corrected potential, the current table-top 
experiment puts a constraint ^ < 14^m [§], which is derived from the results obtained in Ref. In this paper, we 
consider this RS-II model. 

Next, we explain a non-trivial constraint on £ in the RS-II model that can be derived from the astrophysical obser- 
vation of BHs. So far, theoretically no brane-localized BHs larger than £ have been constructed either analytically or 
numerically. By employing the AdS/CFT correspondence [111, [13] to the brane-localized BHs, it has been conjectured 
that such BHs cannot be static [H, (see also Ref. as an application of the AdS/CFT correspondence to the 
RS-II braneworld scenario). Applying the AdS/CFT correspondence, a five dimensional BH is considered to be dual 
to a four dimensional BH associated with CFT fields. The latter system should evolve via Hawking emission from 
the BH. According to the dictionary of the AdS/CFT correspondence, the number of degrees of freedom of CFT is 
as large as 

9. = l5/^-l0''(-^y . (3) 



G \10fim 

This enormously lar ge f actor enhances the BH evaporation rate considerably. For a BH with mass Af , the evaporation 
rate is evaluated as [l6| 

where Cj^j was defined as the coefficient of the mass loss rate for later use. This leads to the estimate of the lifetime 
of the BH, 
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yr. (5) 



Emparan et al. [Ta | pointed out that primordial BHs, if detected, may put a strong bound on £. Several constraints 
have already been obtained from astrophysical BHs by (i) estimating their masses and ages and by (ii) measuring the 
orbital decay rates of BH binaries. For the former cases, Psaltis [131 estimated the lower limit on the age of the BH 
in the X-ray binary XTE J1118-I-480 to be llMyr, which leads to ^ < 80/xm. Gnedin et al. [l^ also estimated the age 
of the BH in the extra-galactic globular cluster RZ2109, and obtained a conservative bound £ < lO^m. For the latter 
cases, Johannsen et al. |l9l. [20j focused on the upper boimd on the orbital decay rates of the X-ray binaries A0620-00 
and XTE J1118+480, and placed bounds £ < 161/^m and £ < 970/im, respectively. The inclination i of the binary 
A0620-00 was estimated as i 51.0° ± 0.9 by Cantrcll et al.~^^d it leads to the BH mass M = 6.6 ± O.25M0. 
Since this is almost 2 times smaller than the one assumed in Refs. |19l.[20|. this new results may put a slightly stronger 
constraint. Recently, Simonetti et al. found that if the orbital decay rate of a BH-pulsar binary is detected in future 
with the same accuracy as 30 years observation of PSR B1913-I-16 the 5 a upper bound on £ becomes 0.17/xm [23j . 

It is also possible to constrain £ from future gravitational wave observations. Inoue and Tanaka [l^l derived 
the leading correction to the gravitational wave (GW) phase due to the modification in the gravitational potential 
mentioned in Eq. By detecting GW signals from sub- lunar mass BH binaries with the 3rd generation GW 

interferometers, they obtained a rather weak upper bound on £. 

Recently, McWilliams [l[ estimated the possible constraints on £ in future by using the Laser Interferometer Space 
Antenna (LISA) [2^ from 2 different observables; (i) the event rate of extreme mass ratio inspiral (EMRI) and (ii) 
GW signal from a galactic BH/neutron star (NS) binary. From the former estimate, it was claimed that £ can be 
constrained as £ < 6/zm for a (5 4- 10^)Mq binary if the predicted EMRI event rate in General Relativity (GR) [2^ is 
correct. From the latter estimate, a monochromatic binary signal of a galactic {2+5)Mq BH/NS binary at / = lO^'^Hz 
puts a constraint £ < 22/j,m (see Ref. [131 the relativistic stars in RS-II model), assuming that it is in the inspiral 
phase. However, in order to probe the mass loss effect more robustly, we need to determine the change in the orbital 
separation, which cannot be made from a monochromatic signal. Namely, we need to detect chirp (or anti-chirp) of 
GW signals from binaries. In detecting these signals from a stellar mass BH/NS binary, the Deci-Hertz Interferometer 
Gravitational Wave Observatory (DECIGO) [l^, [2^ and the Big Bang Observatory (BBO) [s^, both having optimal 
sensitivities at 0.1-lHz, perform better than LISA. 



3 



In this paper, we first derive the correction to the GW phase due to the mass loss effect, whose frequency dependence 



behaves like "-4PN" correction (see Ref. l31[ for a related work). Then, we perform the matched filtering analyses 
and estimate the possible constraints on £ by detecting GWs from BH binaries with LISA and DECIGO/BBO. 
Since high BH/NS event rate of 0(10'^ yr'^) has been predicted for DECIGO/BBO in GR [13, we can obtain a 
stronger constraint by performing statistical analyses. We also estimate the upper bounds on £ from the expected 
detection number of EMRI events for LISA and of BH/NS ones for DECIGO/BBO, extending the previous work by 
McWilliams 

This paper is organized as follows. In Sec.|lTl we review 2 ways of obtaining the constraints on t with LISA estimated 
by McWilhams ^. In Sec. Hill we re-examine the constraints by performing matched filtering analyses. First, we 
derive "-4PN" correction term appearing in the gravitational waveform phase. Next, we describe the basics of the 
Fisher analysis and explain the future planned space-borne interferometers that we use for our analyses. Then, we 
evaluate the constraint obtained from a single event. After that, we perform statistical analyses and show that the 
constraints are improved. In Sec. IIVI we extend the discussions made in Ref. [l| for obtaining constraints from the 
number of detection events. We consider not only EMRI detection numbers with LISA but also BH/NS ones with 
DECIGO/BBO. Finally in Sec.|Vl we summarize our work and comment on several issues that we did not take into 
account in this paper. We also mention possible future works at the end. We take the present Hubble parameter as 
i^o = 72km s^^ Mpc^^ and the cosmological density parameters as Q,„i — 0.3 and 51a = 0.7. Hereafter, we take the 
unit G = c = \ throughout this paper. 

II. CONSTRAINTS WITH LISA OBTAINED BY MCWILLIAMS 

In this section, we review two ways to constrain the size of extra dimension I using LISA developed by 
McWilliams LISA will detect an almost monochromatic GW signal from a galactic binary composed of a BH and 
a NS. The event rate for such a binary in the LISA frequency range is expected to be 1 yr~^ [3j|. GW emission makes 
the orbital separation a smaller (inspiral) with the orbital decay rate given by (33 | 
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aGW = --= —1 (6) 

O Qj 

where Adt = m + Ad \s the total mass of the binary and /i = mM/Mt is the reduced mass. On the other hand, the 
BH mass loss effect makes a larger (outspiral) at the rate of 

a« = ^^a, (7) 

where the mass loss rate is given by Eq. This is derived from the conservation of the specific orbital angular 
momentum j = \/Mta assuming that the radiation is emitted isotropically in the rest frame of the BH [s^. There 
exists a critical separation Ocrit where acw and hu balance. If the separation is larger than Odit, the mass loss 
effect dominates over the GW emission and the separation gets larger while if a is smaller than Ccriti GW emission 
wins and the separation gets smaller. Typically, a galactic BH binary forms with its orbital period O(days) [1^ 
whose GW frequency being slightly lower than the LISA sensitivity band. Therefore if its signal is detected at 
/ = IQ-'^Hz, it means that GW emission effect is dominating over the mass loss effect at this frequency. The 
inequality a(/ = 10^'^Hz) < Ocrit leads to the constraint 

where the typical BH and NS masses are assumed to be M = hM^ and m = 2Mo, respectively. 

McWilliams also obtained a constraint on I from the average EMRI event rate (72.)EMRr rate is estimated 

as m 



M fbH^V'"^ _3 



Wemri^(y^) (^) Gpc-V-^ (9) 



However, if we take into account the conjectured BH mass loss effect in the RS-II braneworld model, the evaporation 
time (rather than the age of the universe) also affects the event rate. The modified estimate for the average EMRI 
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event rate (7?.)jj becomes 



EMRI 



where the BH hfctime r is given in Eq. ©. In Ref. [H, the author assumed that this event rate obeys the Poisson 
probabihty distribution and estimated a constraint on I that will be obtained if (7?.)emri' the value predicted in GR, 
is actually observed. However, here we point out that it is not the EMRI event rate but the detection number of 
EMRIs that obeys the Poisson probability distribution. We also have to take into account the large uncertainties in 
the event rate estimations. We study these issues later in Sec. IIVI 



III. CONSTRAINTS FROM THE MATCHED FILTERING ANALYSIS 



In Sec. ini we mentioned that the typical galactic BH binary forms with a GW frequency slightly lower than the 
LISA sensitivity band. However, some of them may form with a frequency higher than / = 10~^Hz [36| . Therefore, 
even if we detect a monochromatic binary signal at / = lO^^Hz, we cannot immediately conclude that the binary 
separation is getting smaller. What we need to detect is the changing rate of the separation which leads to the 
variation of the GW frequency. In other words, we need to detect a chirping or anti-chirping signal. In this section, 
we first derive the correction term in the GW phase due to the mass loss effect. Then, we evaluate the possible 
constraint on i from the matched filtering analysis using LISA or DECIGO/BBO. Throughout this paper, we assume 
that the binaries are quasi-circular. We also neglect the spins of BHs and NSs. 



A. Binary Waveforms 



First we study the GW waveform from a binary composed of two BHs with masses M and m (with M > m). From 
Eqs. (HI and ([7]), the orbital separation change due to the mass loss effect becomes 
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Then, from the total separation shift a = acw 
il 96 
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ciH, the GW frequency shift becomes 
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up to 2PN orders. Here = M/^^/a^^^ is the orbital angular velocity of the binary and x 



(12) 

(7rM,/)2/3is 

the squared velocity of the relative motion. We have also introduced 2 mass parameters, the symmetric mass ratio 
T] = ^/Adt and the chirp mass M = Mtrf/^. The first term in the square bracket in Eq. ([T^ represents the leading 
Newtonian term. The second one is the "-4PN" correction term due to the mass loss effect which we take up to 0{f'). 
The rest of the terms represent the higher PN contributions up to 2PN [33| . 
Next, we integrate Eq. Q to yield 
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where the subscript denotes the quantity at the time of coalescence. The expression for m is obtained just by 
replacing AIq with mg. Using this mass formula, we integrate Eq. (1121) to obtain the time before coalescence t(/) and 



the GW phase 0(/) = J^^^'' 27r/di as 
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and 
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respectively. Here, we introduced L = i^/Alf^, xq = {iTAItof )'^^^ and 
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The binaries of our interest satisfy the foUowing conditions, d\nA/dt <JC d(j>/dt and d^(t)/dt^ ^ {dip/dtY , where A 
is the GW amphtude and is the phase in the time domain. Then, using the stationary phase approximation [3^ . 
the gravitational waveform in the Fourier domain is given as 
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In this paper, we only keep the Newtonian quadrupole term for the amplitude and average it over the directions and 
the orientations of the binaries 133, 



A 



1 



M 



5/6 



307r2/3 Dl 



(18) 



where is the luminosity distance. On the other hand, we keep the phase up to 2PN order. This is the so-called 
restricted 2PN waveform. The GW phase in the Fourier space is given as 
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The second term in the bracket is the "-4PN" correction term in GW phase due to the mass loss effect. Notice that 
C in Eq. (|16p gets larger as rj gets smaller. This is because when Mt is fixed, smaller 77 means smaller m, leading to 
(i) larger mass loss effect, and (ii) suppressed GW radiation. Note also that C = when 77 = "'"'^^4^ 0.14. 
For a BH/NS binary, the orbital separation change becomes 



dH - C*,;. 
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and Eq. ([T^ becomes 



(20) 



/ ^ ^7r^/3_A^5/3jll/3 

5 

Then, the coefficient C changes to 



^ (1 - 277) -~ _4 

48 27,3 Mi"" 



(higher PN terms) 
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C = 
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(22) 



B. Fisher Analysis 



The detected signal s{t) ~ h{t) +n{t) contains both the GW signal h{t) and the noise n{t). We apply the matched 
filtering analysis to estimate how accurately we can determine the binary parameters 6 [34 . [38j with future planned 
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space-borne GW interferometers. We assume that the detector noise is stationary and Gaussian. Then, the noise n 
follows the probability distribution given by 



p(n) oc cxp 



where the inner product is defined as 
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and the quantities with tilde are the Fourier components. Snif) the noise spectral density of each interferometer. 
The signal to noise ratio (SNR) p using a detector with Ni^t effective interferometers is defined by [t^] 



p = ^N~^{h]h). (25) 
Given a GW signal s(t), the probability distribution that the parameter 6 is the true parameter set becomes 



p{h{9)\s) oc exp 



where is called the Fisher matrix defined by 
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dh 



Then, the determination error of the parameter 6^ is estimated as 



(r-i). 



int 



(26) 



(27) 



(28) 



Note that (r-i)^^ is the variance of the parameter 9' using a single interferometer when all the other parameters have 
been marginalized. 



C. Detector Noise Spectrum 

In this subsection, we briefly explain 3 future planned space-borne GW interferometers, LISA, DECIGO and BBO, 
and introduce their noise spectra. First, LISA is an all-sky monitor having a quadrupolar antenna pattern psj . It 
consists 3 drag-free spacecrafts in which free-falling mirrors are contained. These spacecrafts form an approximate 
equilateral triangle with the length of each side 5x lO^km. They orbit the Sun 20° behind the Earth with the detector 
plane tilting 60° with respect to the ecliptic. We can arbitrarily choose 2 out of 3 arms to form 1 interferometer. Then, 
we can linearly combine 3 arms to form another interferometer which corresponds to rotating the first interferometer 
by 45°. Therefore this triangular detector contains Mnt = 2 individual interferometers 13911 . 

For the noise spectrum of LISA, we follow the discussion by Barack and Cutler [40|- The non sky-averaged 
instrumental noise spectral density for LISA is given by [Tlj 
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In addition to that, there are confusion noises from the galactic [33[ and the extra-galactic [4l| white dwarf (WD) 
binaries. 



S'gai(/) = 2.1 X 10"^^ 
S^^-s^\f) = 4.2 X 10"''^ 
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FIG. 1: The non sky-averaged noise spectral density for BBO (thick solid curve) and LISA (thin solid curve). We also show 
the amplitudes of GWs from a (10 + W^)Mq BH/BH binary at Dl = 500Mpc (thin dotted line) and a (1.4 + W)Mq BH/NS 
binary at Dl ~ 3Gpc (thick dotted line). Each dot labeled "1 yr" represents the frequency at 1 yr before the binary reaches 
ISCO. 



respectively. Combining all, the total noise spectral density for LISA with a single interferometer becomes 
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, ^rLISA(/) + ^,f' (./) 



(32) 



Here, k ~ 4.5 is the average number of frequency bins that are lost when each galactic binary is fitted out and Tobs 
is the observation period. dN/df is the number density of galaetie WD binaries per unit frequency given by [i^l 



df \mzj 



-11/3 



Hz 



-1 



(33) 



The lower and higher frequency ends of the LISA sensitivity band are taken as /low = 10 ^Hz and /high = IHz, 
respectively. The noise spectrum of LISA is shown as a (blue) thin solid curve in Fig. [TJ together with an EMRI GW 
of (10 + 10^)Mo at Dl = SOOMpc as a (blue) thin dotted line. 

Next, we consider BBO [s^l and DECIGO [1^ [1^. BBO consists of four triangular sets of detectors whose 
configuration is shown in Fig. [2] This corresponds to TVint = 8 individual interferometers. Each triangle has an 
arm- length of 5x10"' km. Its primary goal is to detect the primordial gravitational wave background (PGWB) of 
^GW = 10~^^. Compared to LISA, BBO has an advantage in detec ting the PGWB since the WD/WD confusion 
noise will have a cutoff frequency at around 0.2 Hz [4l| (see Refs. |32l. |43|. |43 | for the discussions of the NS/NS 
confusion noise). In order to detect PGWB, it is necessary to perform correlation analysis [1^, H^. Therefore 2 of 
the 4 triangular detectors are located on the same site forming a star of David. The rest of the 2 detectors are placed 
far apart to enhance the angular resolutions of the source location. DECIGO has almost the same constellation as 
BBO. The main difference is that while DECIGO is a Fabry-Perot type interferometer, BBO is a transponder-type 
interferometer. DECIGO has arm-lengths of 1000km. 

The noise spectrum of BBO is given as follows. The non sky-averaged instrumental noise spectral density for BBO 
is obtained from Ref. 143 as 
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'-'n.BBO 
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(34) 



It has 20/3 times better sensitivity than the one for the sky- averaged sensitivity [13] • As for 5*1^' and S'f^ s'*', we 
multiply them by a factor T = exp{— 2 (//0.05Hz) }, which corresponds to the high frequency cutoff for the white 
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FIG. 2: The configuration of DECIGO and BBO. There are 8 efi'ective interferometers in total, 
dwarf confusion noises. We also have to take into account the confusion noise from NS binaries, which is estimated 

as mm 

(/) ^ 1.3 X 10-4^ f-L) '^'Y — -r^^, — rV^"'. (35) 

where no denotes current merger rate density of NS/NS binaries. Putting all together, the total noise spectral densities 
for BBO with a single interferometer becomes 

+ S:^^-^^\f)T{f) + 0.001 X S^i^if). (36) 

The factor 0.001 in front of 5?^(/) represents our assumption of the fraction of GWs that cannot be removed 
after foreground subtraction [4^|. (With this choice of the cleaning factor, the residual NS/NS foreground noise 
becomes below the instrumental noise.) The lower and higher frequency ends of the BBO sensitivity band are set as 
/low = lO^'^Hz and /high = lOOHz, respectively. The noise spectrum of BBO is shown as a (red) thick solid curve in 
Fig. [H We also show the amplitude of the GW signal from a BH/NS of (1.4 + 10)Mo at Dl = 3Gpc as a (red) thick 
dotted line. DECIGO has been proposed with 3-4 times less sensitive spectrum than BBO. However, this is not the 
fixed design sensitivity and there is a project going on to improve the sensitivity to the same level as BBO. Therefore, 
for the Fisher analyses below, we assume that both DECIGO and BBO have the noise spectral densities shown in 

Eq. dsni)- 

D. Numerical Setups 

In this subsection, we explain how we perform the numerical calculations of the Fisher analyses. We take 

= (lnXo,ln7?o,io,0o,£'L,i) (37) 

as binary parameters, setting tQ = (pQ = for the fiducial values. We evaluate the determination errors of these 
parameters, especially focusing on L, to estimate how strongly one can constrain or how accurately one can measure 
£. We assume that the observation starts Tyr before coalescence. While calculating the Fisher matrix, we perform 
the derivative of h{f) with respect to 9 analytically. We take the integration range of (/;„, /fin) with 

/hi = max{/iow, /Tyr}, /fin = min{/high, /isco}- (38) 

Here fxyr is the frequency at the time Tyr before the binary reaches the innermost stable circular orbit (ISCO), which 
is given as 
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FIG. 3: (left) The upper bounds on i (Eq. (|4T|) by detecting GWs from BH/BH binaries at Dl = 500Mpc with LISA. We set 
mo = IOM0 and vary the values of Mq. The (red) solid, (green) dashed and (blue) dotted curves represent the bounds obtained 
from 1 yr, 3 yr and 5 yr observations, respectively. The horizontal black dotted-dashed line shows the upper bound t < 
obtained from the table-top experiment [y]. SNR for various mass binaries with each observation period is also shown at the 
bottom, (right) We show the 5 yr results for mo = 3 Mq, mo = 5 Mq and mo = 10 Mq with the (red) solid, (green) dashed 
and (blue) dotted curve, respectively. Again, SNR for each binary is shown at the bottom. 



with t{f) given in Eq. p^ . Here, we only take into account the leading contribution from GR. Higher PN contributions 
and the modification from mass loss effect on /ryr are not important for our analysis. 

is the frequency at ISCO. We performed the numerical integration with the Gauss-Legendre routine GAULEG (ioj . 
This quadrature uses the zero points of the n-th Legendre polynomials as the abscissas and the integrand can be 
calculated exactly up to (2n-l)-th order. We take n = 400 for the analyses of single binary GWs in Sec. IIII El and 
n = 2000 for the ones with statistical analyses of 10^ yr""'^ event rate in Sec. IIII Fl 

We use the Gauss- Jordan elimination for inverting the Fisher matrix j49| . In order to make sure that the inversion 
being performed correctly, we first normalized the diagonal components of the Fisher matrix to 1. Then, we take the 
inversion and convert it to the inversion of the original Fisher matrix (see Appendix C in Ref. [50|). We checked that 
our inversion has succeeded by simply multiplying the inversed Fisher matrix with the original one and see how close 
the result is to the identity matrix Sij. 



E. Constraints from Single Binary GWs 

In this subsection, we show the results for the GW signal from a single binary, first when we use LISA and next 
DECIGO/BBO. Here, we take the fiducial value as L = so that we estimate the constraint on £ assuming that GR 
is the correct theory. From Eq. ((28)) . one can easily relate the constraint on £ to the Fisher matrix as 

e<£u^ir-')]!j^Mto, (41) 

with £u denoting the upper bound on £. 
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1. LISA 

Figure [3] shows the constraints on I from unequal mass BH/BH binaries with LISA. We place the binaries at 
Dl = SOOMpc. The detection rate of (10 + 1O^)M0 can be estimated from Eq. © as sal yr~^. First, we set 
Too = IOMq and vary Mq from 10'^ Mq to IO^Mq. The results are shown in the left panels of Fig. [3] for observation 
period of 1 yr, 3 yr and 5 yr, together with the table-top constraint of £ < 14/im Q. We also show the SNR of each 
binary at the bottom. From the upper panel, we see that the upper bounds on ^ from BH/BH binaries with a large 
mass ratio have very weak dependence on Mq, although SNR is increasing as Mq gets larger. This is simply because 
"-4PN" correction term in GW phase is smaller for larger Mto- In fact, we can easily estimate that the upper bound 

— 1/8 

on I from EMRI scales as Mg as follows. 
From Eqs. (HHl-dlll), we get 

^ cx Ar"'C{nM,,f)-^l\nM,f)-^l^ oc M-^/^^-^/V-"/^ (42) 

where we have used C oc ■q^'^ for 770 ^ 1. For the BH/BH binaries considered here, /xyr ranges from 10^'^ — 10^^ Hz. 
In this frequency range, the noise spectral density of LISA can be approximated by a constant. From these estimates 
with Eq. (I?fl) . the L-L component of the Fisher matrix becomes 

« M^^S^'fr^^ (X M;,"%'"'t'^/\ (43) 

where we have used Eq. p9p . If we can neglect the degeneracy between L and the other binary parameters, the upper 
bound on £ follows from Eq. (|4ip as 

C ^ e"™''' = {TLLV'/Hdt^ OC a4^/^%„'^/^^T-i5/i6. (44) 

For the case of EMRI, Mto ^ and fj^- Then, we have 

^(uncor) ^ ^^-l/8^21/16j,_15/16^ (45) 

We find that depends only weakly on Afo- Also, the scaling of y-is/ie jg roughly consistent with the upper left 
panel of Fig. [3] 

Next, we fix the observation period as 5 yr and estimate for too — SAf©, bM^ and IQMq. The results are shown 
in the right panels of Fig. [31 together with SNRs at the bottom. The toq dependence of is consistent with our 
simple estimate above, oc TOo^^^^ ~ TOq '^. From Fig. [31 we see that constraint on t from BH/BH observation with 
LISA is weaker than the one from the current table-top experiment. Also we note that the binaries with relatively 
small Afo may have too small SNRs for applying Fisher analyses [5l|. In this case, thinking of binaries closer to us 
(though the event rate reduces to less than 1 yr~^), we can give a constraint scaling D'. 



^LL oc 



dh 
dL 



r 2 



Snif) 



2. DECIGO/BBO 

Here we examine the estimate of the constraint obtained by using DECIGO/BBO. We fix the binary distance as 
Dl = 3Gpc. In Fig. [H we show the upper bounds on t from equal mass BH/BH binaries with observation periods of 
1 yr, 3 yr and 5 yr. The meaning of each line is the same as in the left panels of Fig. [31 Since we know that smaller 
mass binaries put stronger constraints (because the evaporation times are smaller), we only vary A/o from 3AfQ to 
SOAf©. The results show that, by detecting GWs from a small mass binary, DECIGO/BBO can put a more stringent 
constraint on i than the table-top experiment. For example, 1 yr (5 yr) observation of a (3+3)Af0 BH/BH binary at 
3 Gpc can put a constraint € < lO^m (3.4/im). 

Next, we consider unequal mass BH/BH binaries. First, we fix too ~ lOA/© and vary Mq from 3A/q to 10^ Af© 
with observation periods of 1 yr, 3 yr and 5 yr. The values of obtained are displayed in the upper left panel of 
Fig. [51 Here, the peaks at Afo ^ SOAf© correspond to 770 ~ 0.14 where C ~ 0. For binaries with large Afo, is almost 
independent of Afo, which is the same feature explained in Fig. [31 In the right panels of Fig. [SI we fix the observation 
period as 5 yr and take too = 3Af0, toq = SAf© and too = 10 Af©. Again, only small mass binaries can put a stronger 
constraint than the table-top experiment. 

In Fig. [HI the upper bounds on t from BH/NS binaries with too — 1.4Af0 are shown. We take observation periods 
of 1 yr, 3 yr and 5 yr and vary Afo from 3Af0 to SOAf©. We see that some binaries put stronger constraint compared 
to the table-top experiment and the features of this figure are similar to the ones of Fig. [H 
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FIG. 4: The upper bounds on ^ by detecting GWs from equal mass BH/BH binaries at Dl = 3Gpc with DECIGO/BBO. The 
meaning of each line is the same as the one in the left panels of Fig. [31 Here only 1 curve is shown for the SNR since 1 yr, 3 
yr and 5 yr results are almost indistinguishable. 




FIG. 5: The upper bounds on i by detecting GWs from BH/BH binaries at = SGpc with DECIGO/BBO. The meaning 
for each line is the same as in Fig. O We set mo = lOM© for the left panels and Tobs = 5yr for the right panels. 



F. Statistical Analysis of 10^ events/yr 

It is expected that DECIGO/BBO will detect BH/NS GWs with the event rate of 10^ yr^i. Following Ref. HI], 
we consider the enhancement of the parameter determination accuracy by performing statistical analyses that make 
use of a large number of events. First, we consider a BH/NS binary at a redshift z with the BH mass (at the time of 
coalescence) being Mg. In order to take into account the complete degeneracy between mass and redshift, we replace 
the mass parameters in the Fisher analyses to the redshifted ones mo — >■ rn-oz = (l + 2)mo and Mq — > Mqz = (l-l-z)Mo. 
The way how L is defined in terms of the original total mass, L = jM^^, is unaltered. Then, the total variance a'\i 
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FIG. 6: (Top) The upper bounds on I by detecting GWs from BH/NS binaries at Dl = 3Gpc with DECIGO/BBO. The 
meaning for each curve is the same as the one in the left panels of Fig. O We set the NS mass to mo = IAMq. 



of the parameter £^ is given by double integration of z and Mq as 



^^2 — ^obs 



47r[aor-(z)]^riL(z,Mo) — [crL(^,Afo)A4]-2dzdMo. 
dz 



Here Oq represents the current scale factor and r(z) is the conioving distance to the source given as [32 

dz' 



1 r 

aor{z) = — / 
^0 Jo 



/o ^/nra{l + Z'f + JIa 

T is the proper look back time of the source and ^ is given as [s^l 



dT 



dz i/o(l + z)Vf^™(l + z)3 + r!A 



(46) 



(47) 



(48) 



hL{z,Mo) = /(Mo)/L(Mo)7?.i?(z) [S^l shovifs the BH/NS merger rate at redshift z where TZ is the estimated merger 
rate at today in GR and 



R{z) = 




(49) 



encodes the time evolution of this rate [5 
to satisfy jj^^' 



/(Mo) denotes the mass function of the BH/NS binaries in GR normalized 
j^j^^ f {Ma)dMa ~ 1 where Afo^min and Afg^max represent the minimum and maximum values of Mq in 
the distribution, respectively. /l(A/o) is the reduction rate of the total number of merging BH/NS binaries with 
the BH mass being Mq in RS-H model with L compared to the case of GR. /L(Afo) is estimated by calculating the 
probability that BH/NS binaries merge within the BH evaporation time t (Eq. For the probability distribution 
of the binary merger time tmorg, we use the one shown in Fig. [7] which is a simplified fitting of the result shown in 
Fig. 8 of Ref. ^SGi]. We here assume that this probability distribution does not depend on the BH masses, from which 
we obtain rather conservative results (see Sec. |V|for further discussions). 

We set Too = 1.4M0 and following Fig. 2 of Ref. [s^, we take the fiat mass distribution between Afo,min = SM© 
and Afo,inax = ISAf©. Wc choose TZ = 10^^ Mpc~^yr~^ for our fiducial value (see Sec. |lV]for further details). The 
determination error of £ for the statistical analysis, defined as 



^^(stat) ^ 



(50) 
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FIG. 7: Probability distribution of NS/BH merger time tmerg [3^ . 
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FIG. 8: The detection errors of I (Eq. (|50|) ) with statistical analysis for 1 yr (red solid), 3 yr (green dashed) and 5 yr (blue 
dotted) observations of BH/NS binaries with DECIGO/BBO. The (black) dotted-dashed hue represents A^'^"'*' = i. If A^(=*^*' 
comes below this line, i can be measured. 

is shown against £ in Fig. |5]for various observation periods of 1 yr (red solid), 3 yr (green dashed) and 5 yr (blue 
dotted). As we lower £, A£*^^*^') becomes flat. This is because the BH evaporation time r becomes larger than the 
age of the universe for small £, meaning that the number of BH/NS binaries detected is the same as the one in GR, 
irrespective of the value of I. If IS.i^^^'^^) is below the dottcd-dashcd line that corresponds to A.^^***^'^ = £, £ can be 
detected. From this figure, we see that DECIGO/BBO can measure I down to = 1.8^m (1 yr), I ~ 0.50/im (3 yr) 
and I = 0.33/zm (5 yr). This is about 40 times stronger than the current table-top result Q- On the other hand, if 
the curve of A£(''tat) is above the line of A£("*^') = £, we can not detect such small I but only give constraints on it. 
The constraints on I obtained in this manner are summarized in Table |T] when GR is correct (second row) and when 
RS-II is correct (fourth row). If GR is the correct theory, the upper bound on i becomes the ones shown in the second 
row of Table m Note that the upper bound with 1 yr observation {i =- 1.6/im) is slightly stronger than the detection 
limit oi I = 1.8/im. This difference exists because if RS-II is the correct theory, some of BHs may have evaporation 
time shorter than the age of the universe which reduces the detection rate of binaries, leading to the weaker limit. 
The reason for this difference being small is because Af^''*'^*) is only weakly dependent on the number of binaries as 
^^(stat) ^ (number of binary)^/'''. For 3 yr or 5 yr observation, the upper bound is the same as the detection limit 
value. This is due to the fact that, ii £ = 0.50/xm, the evaporation time of 3M(t,{= Mo,min) BH is about 10^° yr. This 
means that the detection rates are the same in both RS-II and in GR, which leads to the same upper bound. 
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TABLE I: DECIGO/BBO constraints on I from a single (1.4 + 3)A/0 BH/NS binary at Dl =3Gpc (first row) and from 
statistical analyses assuming that GR is the correct theory with (second row) and without (third row) taking eccentricity as 
a parameter. We also show the detection lower limit on I with DECIGO/BBO assuming that RS-II model is the correct 
theory with (fourth row) and without (fifth row) taking eccentricity as a parameter. For the statistical analyses, we assume 
a flat BH mass distribution between 3Mq and 13Mq for BH/NS binaries with the total detection event rate of 10^ yr~^ for 
DECIGO/BBO. If RS-II is the correct theory, we cannot measure I for 1 yr observation with taking eccentricity as a parameter 
since the determination error A^^"*^*' is always larger than the fiducial value of i. 



obs. period (yr) 


1 yr 


3 yr 


5 yr 


single, GR, without ecc. (/im) 


11 


4.3 


3.0 


statistical, GR, without ecc. (^m) 
statistical, GR, with ecc. (/xm) 


1.6 
6.5 


0.50 
2.2 


0.33 
1.4 


statistical, RS-II, without ecc. (/im) 
statistical, RS-II, with ecc. {^.ra) 


1.8 


0.50 
2.7 


0.33 
1.5 
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FIG. 9: The error contour showing 95% confidence level on the le-L plane for a (1.4 + 10)AfQ BH/NS binary at Dl ~ 3Gpc, 
such that the fiducial values lie at the center of the ellipse. It is assumed that the observation period is 1 yr with DECIGO/BBO 
and we take 7^'^'^' = L'*^"*' = for our fiducial values. 

G. Effects of Eccentricities 

Up to now, we assumed that the binaries are quasi-circular. This is because eccentricity e decreases as e cx 
j-19/18 Q gQ ^j^g^^ expected to be negligible by the time when binary GWs come into the observation band. 
However, some binaries, especially EMRIs and intermediate mass ratio inspirals (IMRIs), may have non-negligible 
amount of e even at 1 yr before they reach ISCO. For example, based on the results obtained by Hopman and 
Alexander [5^1, two of the authors of this paper estimated the eccentricity of a (1.4-|-1O^)M0 BH/NS binary typically 
as large as e = 0.026 at / = /lyr (HO]. Therefore we need to include eccentricity into binary parameters for more 
realistic analysis. When averaged over 1 orbital period, the mass loss effect does not change e [2^. Therefore the 
leading contribution of eccentricity to the binary waveform S'i/e is unaltered [50| : 

<5*. = -A(,A^„/)-5/32355,^^-i9/6^ (51) 

where le = (7rMto)^^^^e^/c^^^ is the dimensionless asymptotic eccentricity invariant with /c and Cc being an arbitrary 
reference frequency and the eccentricity at that frequency, respectively. In Fig. [HI we show the 95% confidence level 
contour in I^-L plane for a (1.4-|-1O)M0 NS/BH binary with the fiducial parameters, le = 0, £ = and an observation 
period of 1 yr with DECIGO/BBO. We see that there is strong degeneracy between these 2 parameters. Unfortunately, 
we cannot obtain stronger constraint on £ even if we consider only the right half of Fig. [S) This means that prior 
information of /e > docs not affect the bound on £. In the third row of Table HI we show the constraints on £ 
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FIG. 10: Same as in Fig. |8]but here, eccentricity is included in binary parameters with the fiducial value 1^ ~ 0. 

obtained from statistical analyses assuming that GR is the correct theory and taking eccentricity into account. These 
are about 4 times weaker than the ones without taking eccentricities into account (second row). 

Next, we consider the effect of eccentricity when we try to measure £. We performed the same statistical estimate 
explained in Sec. lIIIFl but this time, we included /g into binary parameters. For the actual computation, we set If, = 
as a fiducial value. The results are shown in Fig. [TOl Since £ and le have strong degeneracy, A^^*'*''*) becomes about 5 
times larger than the ones shown in Fig. |8l In this case, 1 yr observation curve does not cross 

means that when we take into account eccentricity, i cannot be measured with only 1 yr observation period. However, 
3 yr and 5 yr observations will suffice in measuring it with the detection upper bounds of 2.7/im (3 yr) and 1.5 /jm. 
(5 yr), respectively. These are still 1 order of magnitude stronger than the current table-top result. These results are 
also summarized in Table H] for GR (third row) and for RS-II (fifth row) . 

IV. CONSTRAINTS FROM THE EXPECTED NUMBERS OF EMRI AND BH/NS DETECTION 

EVENTS 



In this section, we derive the constraints on £ from the expected numbers of events for EMRI with LISA and 
BH/NS binaries with DECIGO/BBO, extending the discussion in Sec. HI] made by McWilhams [l| with LISA. We do 
not consider BH/BH binaries with DECIGO/BBO here since these merger rates are expected to be about 1 order of 
magnitude smaller than the BH/NS ones [H,!!^]- We define {N) as the averaged number of detection events expected 
in GR and (N)-f ^ as the one in RS-II model. We assume that there is 1 order of magnitude uncertainty in the value 
of (N) (see e.g. ^^)- Due to the mass loss effect, (N)^ is smaller than (N), and we set the ratio rn between these 
two numbers as (N)-^ = th (N) (with rn < 1). Since the detection event number follows the Poisson probability 
distribution (5^, the variance ct^ equals the mean {N)-^. If (N) events are observed and this value (N) is within 5-cr 
of the probability distribution, the following inequality. 



puts the lower bound on th as 



(iV)H + 5aH > {N) , 



2 (N) + 25 - VlOO (A^) +625 



(52) 



(53) 



EMRI Observation with USA 



LISA can observe a (5 -|- 10^)AIq EMRI up to -DL^(max) = 4.3Gpc, assuming that observation starts 1 yr before 
ISCO and SNR threshold is 10 with 2 effective interferometers. Using the EMRI event rate given in Eq. ([9]), the 
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averaged number of EMRI detection events (N) ^mri becomes 

(-^)eMRI ~ ~:r^L, (max) ^obs ('^)eMRI 



3/8 / \ 1/2 



. 3.3 X 10^ f f ^U4^V" f ^ ) . (54) 
V 4.3Gpc J V lyr / V lO^M© J \ m J ^ ^ 

Then Eq. (|53p becomes rn 0.76. This bound is much stronger than the one from the uncertainty in the event rate 
estimation (N)-^-^-^-^, suggesting that the constraints from Poisson statistics are not appropriate in this case. We have 
assumed that the upper bound on within estimation uncertainty wiU be one order larger than the most 

probable value. If 3.3x10^ EMRIs are detected as expected with the theoretical upper bound (^)emri ~ '^■'^ ^ -'^0^, 
the lower bound on rn is set as rn 0.1. From Eq. (fTO)) . this inequality leads to t > lO^yr. Combining this with 
Eq. ([5]), we obtain the upper bound on £ as 

in9 \ 1/2 / \ 3/2 

lO^yr^ ' / \ ' 



When only (TV) = 33 EMRIs are detected with the same theoretical upper bound, the lower bounds on both th and r 
are reduced to 10% of the values mentioned above, leading to € ^ 12/im. These results are summarized in the upper 
half of Table HH 



2. BH/NS Observation with DECIGO/BBO 

For BH/NS event rate with DECIGO/BBO, we do not have a simple analytic estimate like the one in Eq. ([3. 
Therefore we have to rely on the results obtained from population synthesis simulations [IBl- Table 2 of Ref. |36l | 
gives the merger rate of BH/NS binaries as 7?. = 8.1 x lO-'^MWEG" Vr"^ where MWEG stands for the Milky Way 
Equivalent Galaxy. To convert this rate to Mpc~^yr~^. we should multiply a factor 1.1 — 1.6 x 10~^/i72Mpc~'^ [59j . 
This leads to the BH/NS merger rate of 9.0 — 13 x 10~*Mpc~'^yr~^. Following Cutler and Harms [13], we estimate 
the averaged number of detection events (Af)BH/NS '^°i' BH/NS binaries as 

(^)bh/ns === ^ohs / 4:n[aor{z)]'^h{z)—dz (56) 
Jo ""^ 



10^ 



lyr J VlO^'^Mpc^Vr"^ 



where aor{z) and dr/dz are given in Eqs. (|T7|) and (pS]). respectively. This gives the detection rate of (^)bh/ns = 
9.0 — 13 X 10'*, which yields the 5-a bound from the Poisson statistics as 

rn > 0.98 - 0.99. (57) 

Again, this lower bound is much stronger than the one coming from uncertainty in the theoretical prediction for 
(-^)bh/ns Poisson statistics are inappropriate. Therefore it is fair to say that if 9.0-13x10''^ BH/NS binaries are 
detected, we can only put a constraint ru > 0.1 as before. By using Fig. [7l we translate this inequality to the lower 
bound on r as 

r > 6.9 X lOV- (58) 
Using Eq. ^ (or equivalently Eq. ([?5|) '). the upper bound on £ becomes 

r n 1 r,7 \ 1/2 / 71 r \ 3/2 

6.9 X 10 yr\ ( M ^ ' 



<15(— ^) (^) /.m, (59) 



where we assumed BH masses to be M = 5Mq . These results are summarized in the lower half of Table |TI1 We see 
that the constraints obtained in this section are weaker than the ones with statistical analyses in the previous section. 
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TABLE II: The constraints on I from number of detection events. We also show other important parameters used to derive 
the final results. 



Detectors, binaries, masses 


{N) 


r-H 


(yr) 


I 

(^m) 


LISA, EMRI, (5 + 10^)^/0 


3.3 xlO^ 


0.1 


10** 


3.9 




33 


0.01 


10* 


12 


DECIGO/BBO, BH/NS, (1.4 + 5)M0 


9.0-13x10* 


0.1 


6.9x10' 


15 



V. CONCLUSIONS AND DISCUSSIONS 

In thispaper, we obtained the possible upper bounds on the size of extra dimension (. in the RS-II braneworld 
scenario [I by detecting GWs from BH/BH and BH/NS binaries with LISA and DECIGO/BBO. The current table- 
top experiment puts I < 14/im p. It has been conjectured that in the RS-II model, the Hawking radiation is enhanced 
by a factor depending on i [13,llJl- This factor becomes 10^^ for i = lOfim and this mass loss effect may be detected 
by future GW observations [l|. 

First, we derived the "-4PN"-like correction term due to this mass loss effect in the phase of the gravitational 
waveform. It is a negative PN order correction which means this mass loss effect is greater when the separation of a 
binary is larger. Then we performed Fisher analyses and estimated the upper bound on i. The constraints from 
high mass ratio BH/BH observation with LISA are almost independent of the mass of the larger BH Mq, while they 
scale as rrig^^^^ for the mass of the smaller BH mo. We gave some analytical explanations to these behaviors. We found 
that these constraints using LISA are weaker than the current table-top result. On the other hand, DECIGO/BBO 
may put stronger constraint than the table-top experiment by detecting small mass BH/BH and BH/NS binaries. 
For example, when GR is the correct theory, 5 yr observation of a (3-|-3)M0 BH/BH binary leads to ^ < 3.4/im. 
Furthermore, DECIGO/BBO is expected to have very large BH/NS event rate of 0(10^ yr~^)- We found that 
DECIGO/BBO can measure £ down to £ ~ 0.33/^m for 5 yr observation by performing statistical analysis if we know 
eccentricities a priori. This is almost 40 times stronger than the table-top one. When we include eccentricities into 
binary parameters, the constraint reduces to i = 1.4/im for 5 yr observation. When considering actual detection of £, 
the detection limit becomes £ = 1.5/im when including eccentricities. This is still 1 order of magnitude stronger than 
the table-top one. However, table-top experiments are performed model-independently whereas the results obtained 
here can only apply to RS-II braneworld model. Therefore we cannot directly compare these results. Also, since we 
do not take systematic errors on the waveforms and the limitation of the Fisher analysis into account, calculations in 
this paper might be underestimating the bounds (60| . The results are summarized in Table ID 

When performing statistical analyses, we assumed that the distribution of merger time (Fig. ^ docs not depend on 
BH mass but this is not true in reality. The determination error of £ is mainly determined by the number of smaller 
mass binaries. From Fig. 4 of Ref. [36|, we see that smaller mass binaries have relatively smaller binary separations 
at their formations. This is because they obtain large kick velocities when the binary components explode to become 
compact objects so that only those having small separations can survive. Since these smaller separation binaries have 
smaller merger times, it is likely that the peak of the merger time distribution would shift to smaller t,ncrg when binary 
BH mass is smaller. This means that there would be more smaller mass binaries that coalesce within BH evaporation 
times. Therefore our estimates shown in Sec. IIII Fl are conservative in this sense. 

Next, we derived upper bounds on £ from the analyses of expected detection event numbers. We extended the 
discussions in Ref. including large uncertainties in the event rate estimations. We found that if the most likely 
number of EMRIs in GR is detected, LISA can put £ < 3.9/Ltm. BH/NS binary observations with DECIGO/BBO give 
1 order of magnitude weaker constraints. These results are summarized in Table HIl 

Unfortunately, strong constraint on £ cannot be obtained from ground-based interferometers. This is because 
the optimal frequency range is higher compared to the space-borne ones, meaning that "-4PN" correction has less 
contribution. To give a concrete example, we estimated £u from the observation of a (3+3)Mp) BH/BH binary at 1 
Gpc with the third generation ground-based interferometer named Einstein Telescope (ET) |6l| . Its noise spectrum 
is given in Ref. [g^]. We found that £u = 5.6 x lO'^fim with SNR=93, which is much weaker than the one from the 
table-top experiment. Constraints are even weaker for the second generation detectors such as advanced LIGO [63j . 
advanced VIRGO [6| and LCGT 

In this paper, we focused on the correction coming from the mass loss effect. However, there is another one 
appearing in the gravitational potential (see Eq. In Ref. [gjj, Inoue and Tanaka obtained the 2PN correction 

term to the phase of gravitational waveform due to this potential correction. We added this 2PN term and performed 
Fisher analyses but found that our results do not change. 
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In this paper, we performed sky-averaged analyses for simplicity. To take the directions and orientations of the 
binaries into account, we need to perform Monte Carlo simulations similar to the ones performed in Refs. j37ll48l. [sof . 
Also, we neglected the effects of BH spins. The leading contribution of spins appears at 1.5PN with an extra parameter 
/3, the spin-orbit coupling. When we add {3 to the binary parameters and perform Fisher analysis for a (1.4-t-lO)Af0 
BH/NS binary at Dl = 3Gpc with 1 yr DECIGO/BBO observation, we obtain £ < 61/im. Compared to the one 
without including /? (shown in Table |T] as £ < 51/im), the constraint weakened only slightly. Since the mass loss 
correction takes -4PN frequency dependence and spin correction is 1.5PN order, the degeneracy between these two 
is weak. Furthermore, when we include the effect of precession, we expect this degeneracy to be solved [66| so that 
the inclusion of spin parameter would not affect our results much. The Hawking radiation also changes the BH spins. 
When taking BH spins into account, we may have to consider the spin down effect due to the Hawking radiation. 

In reality, the change in the BH mass M is caused not only by the enhanced Hawking radiation but also by the 
accretion onto the BH. We need to clarify whether these two effects can be distinguished by Fisher analysis. The 
mass loss rate due to the former effect can be expressed from Eq. ([4]) as 



2 / „ \ 2 



= -0.25M.,,^^j , (60) 

where M^dd is the Eddington accretion rate written as 

Medd = 2.2 X 10-« (y^) ^^ey^'"'- (6^) 

This shows that if i is small, the mass loss effect might be masked by the effect of mass accretion onto the BH if 
the latter is near Eddington accretion rate. However, the binaries that suffer from accretion near this rate will be 
rare. Even if the accretion rate is as high as this rate, in general matter accretion goes through accretion disk. In 
this case the conservation of specific angular momentum does not hold and Eq. (|lip gets modified, leading to the 
contribution differing from "-4PN". Therefore we expect that these two effects can be distinguished. Also, since M 
dependence of the mass loss rate is different between these two effects, it may be possible to separate them. For 
example, quasi-stationary accretion rate Mace is proportional to the area of the BH yielding Afacc fx [l^, l67l - [69j 
and this has different M dependence compared to Eq. Q. Furthermore, since the properties of BH accretions are 
intrinsic and not universal, we should be able to reduce these effects statistically. These issues are left for future 
works. 
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